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Abstract 

We study potential signals of neutralino dark matter indirect detection by neutrino telescopes in a wide 
range of CMSSM parameters. We also compare with direct detection potential signals taking into account 
in both cases present and future experiment sensitivities. Only models with neutralino annihilation into 
gauge bosons can satisfy cosmological constraints and current neutrino indirect detection sensitivities. For 
both direct and indirect detection, only next generation experiments will be able to really test this kind of 
models. 



^ ! 1 Introduction 
O 

""-■^ ' Our present understanding of the universe is described in the framework of general relativity and cosmology. 

I The densities of its components are related by jjj]: 



where a is the scale factor. This equation and current experimental results suggest and focus on a fiat 
universe with the following parameters 12] : 



Cosmological constant: Qa = 0.7 ± 0.1 

Matter: Q. mat = 0.3 ± 0.1 

baryonic matter: Qb = 0.04 ± 0.01; including fl v i s < 0.01 

cold dark matter: Qdm = 0.26 ±0.1 

Relativistic components: 0.01 < fl re i < 0.05 

neutrinos: 0.01 < fi„ < 0.05 

photons: fi 7 = 4.8± ;| x 10" 5 

Hubble's constant: h = H /WQ km' 1 s" 1 Mpc -1 = 0.72 ± 0.08 

In the framework of the Minimal Supersymmetric Standard Model (MSSM) §J |, |,||, @, the lightest 
supersymmetric particle (LSP) is typically the lightest of the neutralinos xi (= x)> X2, X3> X4> the mass eigen- 
states of the neutral gauge and Higgs boson superpartners (B, W 3 , H^, H^). In the rest of this paper, the 
LSP is thus assumed to be the lightest neutralino xi and is called more generically \- Assuming R-parity 
conservation (R = (— l) s+i+2S ), the neutralino is a good stable candidate for cold dark matter m. In this 
context, all sparticles produced after the big-bang give a neutralino in their decay chain, leading to a relic 
bath of x i n the present universe. These neutralinos could be observed via direct detection (x interaction 
with a nucleus of a detector), or indirect detection of their annihilation products (u, 7, p, D, e + ). We will 
mainly focus in this paper on v indirect detection. The relic neutralinos are gravitationally captured by 
massive astrophysical bodies and accumulated at the centre of these objects by successive elastic scatterings 
on their nuclei. The captured neutralino population annihilates and gives rise to neutrino fluxes which could 
be detectable in neutrino telescopes like Amanda/Icecube, Antares, Baikal. 
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2 Neutralino in the CMSSM 



In the MSSM, the mass matrix of neutralinos in the basis (—iB, —iW 3 , H^, HQ) 



I Mi — mz cos j3 sin 9w mz sin /3 sin 9w \ 

M2 mz cos /3 cos 8\y —mz sin /3 cos Qw 

—mz cos fismOw mz cos /3 cos #w — ju 

V mz sin f3 sin #w — fiz sin /3 cos 6w —fi / 



(2) 



where Mi, M2 are mass term of (7(1) and 5(7(2) gaugino fields, /i is the higgsino "mass" parameter and 
tan/3 =< H u > / < //<j > is the ratio of the neutral Higgs vacuum expection values. The neutralino 
composition is: 

x = N 1 b + N 2 W 3 + N 3 H° d +N i Hl (3) 

and we define its gaugino fraction as gf ra c = \Ni 2 + I-/V2I 2 . 

In this model, the introduction of soft terms in the Lagrangian breaks explicitly supersymmetry, leading 
to a low energy effective theory with 106 parameters. The MSSM is therefore not very predictive, and a 
non biased exploration of its parameter space is not possible. As a first step, we will therefore as usual 
concentrate on gravity- mediated supersymmetry breaking in supergravity |}| inspired models, with Grand 
Unification of soft terms at Egut ~ 2.10 16 GeV parameterized by mo (common scalar mass), mi/ 2 (common 
gaugino mass) and Aq (common trilinear term). To geth er with tan (5 and sgn(jj,), these define a 5 parameters 
constraint MSSM (CMSSM) or mSugra model @,]ll], [I§|, from which the 106 parameters can be deduced 
through renormalization group equations (RGE). 

Due to the large top Yukawa coupling, renormalization group evolution can drive m 2 H J\Q EWSB and/or 
m ^i d \QEwsB to negative values, so that the electroweak symmetry breaking (EWSB) originates purely in 
quantum corrections, which realizes radiative electroweak symmetry breaking. Minimization of the scalar 
potential at the electroweak breaking scale Qewsb yields the condition: 

1 2 m H d \QEWSB ~ m HjQ E WSB tan2 P 2 1 2 1 2 1 /,n 
o m z = 7 T7k \ f 1 \Qewsb ~ ^h u \Qewsb M \Qewsb j W 

2 tan p — 1 tan/3>5 

where usually Qewsb is taken as y/m^ m,- 2 ^| to minimize one loop corrections. 

Such mSugra models offer the theoretical advantage over generic MSSM models that problems such as 
Landau poles, charge and color breaking (CCB) minima are partially addressed when dealing with RGE. 

In mSugra, the lightest neutralino can exhibit two different natures, depending on the input parameter 
values O, ^5) : 

- an almost pure bino-like neutralino for low mo, as the RGE drive M\\q ewsb ~ QA1M\\gut = 
0.41mi /2 << \h\q ewsb and M 2 \q bwsb ^ 0.83Mi|gut = 0.83mi /2 << \^\q ewsb \ 

- for mo > 1000 GeV, the neutralino picks up some higgsino mixing for tan /3 > 5 as the increase of mo 
drives m 2 Hu to less negative values so that both \m% u \ and (via eq. ^) decrease. One can then have 
\iAQewsb ~ Mi\q ewsb depending on mi/ 2 - When \fi\ is too small, EWSB cannot be achieved. 

In addition, the very low mo values can lead to models with tachyonic sfermions (low my 2 ) or a slepton 
LSP. These behaviors are shown on figure ^in the (too, w.1/2) plane. The neutralino mass is ~ 0.4mx/2 but 
is affected by low \i values at high mo. Unless explicitly quoted we use m t = 174.3 GeV. 

In this paper, the RGE, SUSY spectrum and potential minimization are derived using the SUSPECT 
program fig, |l7| version 2.005 which includes 2 loops RGE the tadpole method to minimize the scalar 
potential Jl9| | and all radiative corrections in the mass spectrum j2o| . The neutralino relic density, including 
XX + an d XX2 coannihilations, and detection signals are calculated by the DarkSUSY (DS) package J2I], ^ 
which correctly implements s-channel poles and thresholds p3[| . Our interface (available from nezri@in2p3.fr, 
submitted for inclusion in DS) shortcuts the DS spectrum recalculation, to comply with the SUSPECT 
spectra which have been checked E3] against SOFTSUSY & 



3 Branching ratios in neutralino annihilation 

The leading processes in neutralino annihilation at rest are shown in figure ^, where the last two processes are 
proportional to the higgsino/wino fraction of \ an d x! + '- I n mSugra models, one can typically distinguish 
different regions of annihilation branching ratios in the (too,TOl/2) plane, for almost any values of tan/3(> 5) 
(see figure fa). 

1 These calculations use a neutralino local density p x = 0.3 GeV/cm 3 and a Maxwellian velocity distribution with vo = 
220 km s" 1 
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A0=0 ; tan(P)=45 ; \i >0 
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Figure 1: a) Gaugino fraction, b) value of fx and c) neutralino mass in the (mo,m 1( / 2 ) plane. Region where \ is not the 
LSP or EWSB is not achieved are also indicated. 
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Figure 2: Leading channels in neutralino annihilation at rest, with the parametric dependance of their amplitude, Ni 
being defined in equation ^ M . 
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Figure 3: Typical regions of neutralino annihilation. 



Above top thresholds: i.e. m x > m t0 p <^ wii/a ~ 400 GeV or slightly higher in the mixed higgsino 
region. 

Region I: for low mo (< 500 — 1000 GeV), the neutralino is almost purely bino-like. Despite a tiny 
higgsino fraction, its annihilation via the Z exchange amplitude, being proportional to the fermion mass, is 
surprisingly important in the top channel \X tt cc (mtm x /m%) 2 N^, 4 y Another potentially important 
annihilation process is \X ~* bb, whose amplitude \A^\ cc ta ^ m t ( I^jl ) 2 jVi N 3 ^ is strongly enhanced for 
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high tan/3, both from the explicit dependence, and because m,A decreases when tan/3 increases. In the 
region of interest, their ratio \Af^/A^\ accordingly drops from ~ 20 — 50 for tan/3 = 10 down to ~ 2 — 5 for 
tan /3 = 45. Even taking the smaller ti phase space (a factor of 0.7) and the other subleading bb processes into 
account, the bb domination thus calls for an explanation on figure ^ and even more on figure [13| The only 

way to account for this domination is to recall the opposite sign of the XX ~~~ > tt an d XX ~* ti amplitudes, 
thus allowing for a cancellation provided their order of magnitude match, which we checked analytically. 
The opposite happens for the annihilation into bb: the Z, A and sbottom exchange all add up, and from 
the behaviour of the non-negligible amplitude, we expect in the region I an important total annihilation 
cross section 0V_ X f° r m S n values of tan /3. 

Region II: when mo increases (< 1500 — 2000 GeV), tua increases and, for intermediate tan/3 values, the 

A — t — 

XX — * bb amplitude decreases. Since sfermion masses also increase, the XX — * tt channel no longer cancels 

Z — A ~ 

XX — * tt, which becomes the dominant process. However, for high tan/3, XX — * bb is enhanced and stays 
dominant. Globally, since all scalar masses have been increased by mo, and the higgsino fraction is still 
small, the neutralino annihilation cross section cr x _ x is smaller than in region I. 

Region III: increasing further mo (> 2000 — 2500 GeV), for any value of tan /3, one finally approaches the 
mixed higgsino-bino region. The XX ~* it channel amplitude Af^ oc ( "^T* )-^3(4) then dominates all other 
processes, which are suppressed by the increase of scalar masses. So we are left with a ti region parallel to 

A ' 

the highest higgsino fraction isocurves. The mo value separating this region from the previous XX — * bb 
region depends on tan (3. All in all, the neutralino annihilation cross section o" x _ x is strongly enhanced by 
Z exchange near the no-EWSB boundary. 



A0=0 ;tan(P)=45 ; n>0 
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Figure 4: Dominant branching ratios of the neutralino annihilation in the (mo,m 1 /2) plane 



Between top and W/Z thresholds: 
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Region IV: since m x < m top , the neutralino annihilates almost only into bb. Even if some \X — * bb 

A — 

occurs, the XX — * bb amplitude is dominant. 

Region V: increasing mo (and m^), XX ~* bb remains dominant and this process is still quite efficient 

A — Z — 

for high values of tan /3. For intermediate values of tan /3, XX — * bb and XX — * bb both occur but their 
amplitudes are small and the total annihilation o" x _ x is not efficient. 

Region VI: increasing further mo disfavours scalar exchange, but even small, the higgsino fraction allows 

x + 

XX — W + W~ and XX ^* to dominate and enhance the total annihilation cross section cr x _ x . Again, 
the mo values delimiting the boundary with region V depend on tan/3 (via mj). 
Below W/Z thresholds: 

Region VII: the dominant process is XX ~ * bb via A and/or Z exchange, depending on mo, tan/3 and the 
higgsino fraction. 

This analysis is illustrated on figures and Fiji showing the four most important branching ratios for 

A A - 

tan/3 = 45 and 10 (the omitted process XX — * TT behaves as XX — * bb but with a smaller amplitude due 
to the m T /mb ratio). It further offers a qualitative understanding of the relic density picture (figure ^J). 
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Figure 5: Neutralino relic density in the (mo,tni/2) plane. Grey area indicates the excluded models by current acceler- 
ators constraints as discussed in the text. 




The annihilation is efficient for low values of mo (depending on tan /3) and for a mixed neutralino along the 
no-EWSB boundary. This gives the "V" (or "U" ) shape of the relic density profile for large (or small) tan j3. 
According to the current cosmological parameter values, we take the neutralino as an interesting cold dark 
matter candidate if 0.025 < Q,h 2 < 0.3. Figure [s] also shows the region in the (mo,™^) plane excluded 
by the experimental constraints from the Particle Data Group 2000 Pqj implemented in DarkSusy, that we 
have updated with: 

m x+ > 104 GeV; m f - > 100 GeV for / = ti, bi, I*, v, m- a > 300 GeV; m<j 12 > 260 GeV for q = it, d, 5, c and 

-6 x lO" 10 < a„(SUSY) < 58 x 10~ 10 J§7|, ||. 

As noted in previous studies Jis], the muon anomalous moment a^and b — > branching ratio 
constraints strongly favour fi > 0, to which we restrict in what follows. In figure |j| the grey tail at large mo 
is excluded by the chargino bound; it directly bites into the region relevant for neutrino indirect detection. 
Less relevant is the exclusion of small mo and m\n values, which comes both from the Higgs mass limit and 
the b — *■ S7 branching ratio, as calculated by default in DS. The range BR(b -> si) = 1 -> 4 x 10~ 4 chosen 
by default in Darksusy may seem too low in view of latest CLEO and Belle measurements ]2^] . However the 
leading order calculation [B0J implemented underestimates the SM value to 2.4 x 10 -4 , while next to leading 
order corrections give 3.6 x 10~ 4 J3l| so that this range should roughly correspond to 2.2 — > 5.2 x 10 -4 , 
excluding a bit more than the range 2 — > 5 x 10~ 4 chosen for instance in pM- We have checked that replacing 
the implemented Higgs limit J32| by an aggressive version of the latest limit J33| (mj > 114 GeV for all 
sin (/3 — a)) only excludes a few more points around (mo = 1000, mi/j = 150) where neutrino fluxes are 
beyond reach. 

In mSugra, XX + an d XX% coannihilations (included in DS) happen only in the mixed neutralino region, 
decreasing further the relic density, yf coannihilation happens for low values of tan /3, for which there is 
no mixed region, xt coannihilation J34|, [n], ^] happens for high values of Ao. In both cases, sfermion 
coannihilations (missing in DS) are relevant to lower the relic density in regions of large m x , which as we 
will see, are beyond reach of indirect detection. We therefore do not expect their proper inclusion to change 
our conclusions. 
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4 Neutralino-proton cross section: capture rate and direct 
detection signals 

Capture: 

If present in the halo, relic neutralinos must accumulate in astrophysical bodies (of mass Mb) like the Sun 
or the Earth. The capture rate C depends on the neutralino-quark elastic cross section: cr x _ 9 . Neutralinos 
being Majorana particles, their vectorial interaction vanishes and the allowed interactions are scalar (via 

\q ' > XI m * channel and xi XI m s channel) and axial (via xi ~* XI m t channel and XQ XI m s 
channel). Depending on the spin content of the nuclei N present in the body, scalar and/or axial interactions 
are involved. Roughly, C ~ ^ 12n MbfN m CT ^ — < Vesc >n F(v x ,v esc ,m x ,m,N), where p x , v x are the local 
neutralino density and velocity, /jv is the density of nucleus N in the body, ctat the nucleus-neutralino 
elastic cross section, w esc the escape velocity and F a suppression factor depending on masses and velocity 

2 

mismatching. The neutralino capture is maximized for m x ~ mjv as — — — ~ — — — = m x m N . g 
much more efficient in the Sun than in the Earth as Mq >> Mgj . 

For the Earth, scalar interactions dominate. By increasing mo from its low value region, sfermions and 
H exchanges first decrease, and the cross-section rises again when approaching the mixed higgsino region 
(see figure ^a). The capture rate is resonant for m x ~ 56 GeV around the iron mass. 

For the Sun, the spin of hydrogen allows for axial interaction, which are stronger due to the Z coupling. 
The latter depends strongly on the neutralino higgsino fraction and is independent of tan f3, so the cross- 
section follows the higgsino fraction isocurves as can be seen by comparing figure pp and 0. 




To summarize, due to the large solar mass and the spin dependent neutralino-quark cross section, the 
storage of neutralinos is much more efficient in the Sun than in the Earth. 
Direct detection: 

The elastic cross section of a neutralino on a nucleus also depends on o x - q , the nucleus mass number A and 
its spin content. Depending on the chosen nuclei target, current and future direct detection experiments are 
sensitive to the scalar coupling (CDMS (Ge) ||], Edelweiss (Ge) |||, Zeplin (Xe) ||]) or to the axial 
coupling (MACHe3 ( 3 He) [0). Comparison between experiment sensitivities and a set of mSugra models 
are shown on figure bk. 



5 Neutrino indirect detection 

As XX ~ * v v i s strongly suppressed by the tiny neutrino mass, neutrino fluxes come from decays of primary 
annihilation products, with a mean energy E v ~ to . The most energetic spectra, called "hard" come 
from neutralino annihilations into WW, ZZ and the less energetic, "soft", ones comes from bb [El]]. Muon 
neutrinos give rise to muon fluxes by charged-current interactions in the Earth. As both the v charged- 
current cross section on nuclei and the produced muon range are proportional to E v , high energy neutrinos 
are easier to detect. Considering that the population of captured neutralinos has a velocity below the escape 
velocity, and therefore neglecting evaporation, the number N x of neutralinos in the centre of a massive 
astrophysical object depends on the balance between capture and annihilation rates: N x = C — CaN x , 
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where Ca is the total annihilation cross section cr^_ x times relative velocity per volume. The annihilation 
rate at a given time t is then: 



Va = ^C A Nl = - tanh 2 ^CClt (5) 

with Fa ~ % = cs te when the neutralino population has reached equilibrium, and Fa ~ \C 2 CAt 2 in the 
initial collection period (relevant in the Earth). So, when accretion is efficient, the annihilation rate does 
not depend on annihilation processes, but follows the capture rate C and thus the neutralino-quark elastic 
cross section. 




For the Sun, neutralinos do not reach complete equilibrium in the whole (mo, Wu/2) plane studied, which 
contrasts with (29) based on Neutdriver j^] and different RGE, probably yielding lighter (pseudo-)scalars and 
lower \i. We find that the annihilation is on average less efficient (as attested by our smaller region with an 
acceptable relic density), which results in a filling fraction \j CCa^q smaller by a factor of fO to fOO. This 
has little effect on u(fi) fluxes for high mo or small m 1 / 2 < 600 GeV values, as equilibrium is nevertheless 
reached and fluxes follow essentially the higgsino fraction and the spin dependent <7 X ^!™ isocurves of figures [l]a 
and|^. For low mo however, the equilibrium fluxes would drop with mo and C, but the smaller capture rate 
hinders equilibrium (e.g. ^JCCa^q ~ 1 for mo = 500, mi« = 800) and Ya feels the increasing annihilation 
cross-section. This effect is stronger for high m.1/2 values where the neutralinos are heavier and more bino- 
like, and where capture is smaller. The effect is also larger when tan/3 (and thus 0" x _ x ) ^ s l° w ' as incomplete 
equilibrium makes f(n) fluxes sensitive to cr^_ x , making the W and top thresholds more conspicuous. 

A0=0 ; tan( P)=45 ; fi>0 
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Figure 8: v fluxes from the Earth (a), the corresponding /i fluxes _E M > 5 GeV theshold (b) and their ratio (c). 
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For the Earth, neutralinos are not in equilibrium. Neutrino fluxes depend both on C 2 and annihilation, 
giving an enhancement in the low and high mo regions where fluxes are boosted by annihilation (see figure ^ 

< o~l P 2v, the capture rate and v fluxes from the Earth are much smaller 



andH). Since Mq < Mq and o x 
than from the Sun. 



seal 
V 



Comparing v fluxes and /x fluxes, the v — * fj, conversion factor increases with m 1 / 2 due to the increase of 
m x leading to more energetic neutrinos. This ratio also follows the annihilation final state regions described 
in section 3 above. Indeed, spectra from WW, ZZ and to a lesser extent ti are more energetic than bb 
spectra, so neutrino conversion into muon is more efficient in the mixed bino-higgsino region above W and 
top thresholds (see figure |^ and ^) . 
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Figure 9: a) Direct detection experiment sensitivities and a set of models in the (o x 



plane. The models excluded 



by current accelerators constraints are not displayed, b) Indirect detection experiment sensitivities in the (fi fluxQ , m x ) 
plane with — 5 GeV threshold of the same sample as a). Dotted, dash-dotted and dashed curves are respectively the 
Macro, Baksan and Super-Kamiokande upper limits. Solid lines indicate the future Antares and Icecube sensitivities. 



Comparison with indirect detection experiment sensitivities: Muon fluxes coming from neutralino anni- 
hilation in the Sun are shown in figure Bp. MSugra models predicting a good relic densit y o f neutralinos 
give neutrino/muon fluxes which can be as high as the current experimental limits (Baksan |4l|| , Macro 0], 
Super Kamiokande ^JJ). The 0.1 km 2 Antares detector will explore further the interesting parameter space 
§44]. Next generation neutrinos telescopes (Icecube, Antares km 3 ) will be much more efficient to test such 
models, eg Icecube expected sensitivity ~ 10 2 /i km 2 yr _1 from the Sun ]4q] . 

Indirect vs Direct detection: a high neutralino-proton cross section is efficient in both direct and indirect 
detection (via capture). The large mo mixed higgsino-bino region pointed out in j^, |i[ favours both direct 
and indirect detection due to the enhancement of a x - p and o~ x _ x . Both enter in the indirect detection which 
is moreover favoured by the production of more energetic neutrinos in WW, ZZ and ti decays, leading to 
better conversion into muons (figure ^). This mixed region, which has a good relic density, is very attractive 
for neutrino indirect detection signal. 

Hovewer, no current experiment is able to really test such models. Figure [ll] shows the region of the 
{mo,mi/2) plane which can be explored by future direct detection experiments and neutrino telescopes. 



6 Variations 

In this section, we discuss the robutness of the indirect detection picture described above, with respect to 
variations in Suspect 2.002 of 1) mSugra input parameters like Ao or tan/3, 2) experimental uncertainties on 
the top quark mass and 3) theoretical uncertainties on the scale of electroweak symmetry breaking Qewsb 
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A0=0 ; tan(P)=45 ;n>0 




Ao: Varying Ao away from does not change too much the analysis above. Annihilation branching ratio 
regions, as well as the no-EWSB boundary and slope can be slightly displaced. For moderate tan /3 ~ 20, 
quite large and negative values of Ao (eg —800 GeV) can enhance r + r~ annihilation channel along the f LSP 
region (i.e. low mo) due to splitting in the f mass matrix. For such values, the fx coannihilation [^(| could 
also affect the relic density, but in this region, neutrino fluxes are too small to be detected anyway (see figure 
[[§). The same happens with the tt channel for 400 < m 1/2 < 600 and m < 700 (~ \A \). The t X Q 0, fj| 
coannihilation can also happen for Ao ~ —2000 GeV, but the usual cosmologically interesting region is re- 
duced (due to the negative squared mass arising in sfermion matrices at low mo and because the mixed 
region is pushed to higher values of mo). In addition, this coannihilation region produces neutrino fluxes 
from the Sun around 10 _1 /x km -2 yr _1 , too low to be detected. In view of the latest results on cosmological 
parameters, figure |l2]a shows models with 0.03 < Qh 2 < 0.3 in the (mo, m^) plane for different values of Ao, 
as well as models giving fj, fluxes from the Sun larger than 10 2 (b) and fO 3 fj, km -2 yr _1 (c). Latter models 

are confined in the mixed higgsino-bino region. In addition, only those models with mw < rrix < m « gi vm S a 

x + 

hard spectrum of neutrinos via %X — ^ W + W~ and XX — ^ ZZ can simultaneously satisfy 0.1 < flh 2 < 0.3, 
while yielding a high muon flux. 

tan /3: the branching ratio picture does not qualitatively change for different values of tan j3. Low values 
of tan/3 can add a \X ~~* ti region for low mo and open a f\ coannihilation region along the f LSP 
region, but \x fluxes are small and such values of tan f3 reduce the mixed neutralino region which is the most 
interesting for indirect detection. As discussed previously, the boundaries of branching ratio regions move 
with tan/3 (compare figure [l3] and ^J). When tan/3 grows, the cosmologically motivated region is wider due 
to the annihilation enhancement via A exchange (in regions I /II and IV/V). Applying the conservative cut 
0.03 < Q,h 2 < 0.3, all values of tan/3 provide interesting models for neutrino indirect detection (see figure 
[li]). For too large values of tan 0, the bb region is extended and leads to soft neutrino spectra. In this case it 
becomes more difficult to conciliate a tighter limit on the relic density with current experimental sensitivities. 
Values of tan /3 ~ 10 — * 35 give more models satisfying both relic density 0.1 < Q,h 2 < 0.3 and high ^//ifluxes. 

EWSB scale: Qewsb is the scale at which the running of the soft SUSY breaking terms is frozen and 
the minimization of the scalar potential is achieved. When the full potential is minimized within a multi- 
scale RG treatment, Qewsb should not affect physics. Hovewer, the one-loop single scale effective potential 
does depend on Qewsb, and the tree- level relations (equation ft) only hold in terms of running parameters 
at Qewsb', if Qewsb is chosen ~ i/ m i 1 m t 2 [j0§> ^ ne potential can be minimized pertubatively. Varying 
Qewsb away from this value is therefore not physical, but gives a hint on the theoretical uncertainties 
associated with a single scale RG and potential minimization approach. When Qewsb is lowered, the RG 
running drives /i to smaller values. This leads to a wider region of the (mo,mi/2) plane where EWSB 
cannot be achieved. In the remaining allowed region, the lower values of /i slightly increases the higgsino 
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fraction of the neutralino and the annihilation cross-section, giving more cosmologically acceptable models. 
Figure |l5| and Ft] show the relic density and /i fluxes from the Sun and the Earth in the (mo, m 1 / 2 ) plane 
for Qewsb = 2 \/ m 'i m * i ant ^ | V m *i m *i ' ^ or Q ew sb — gV m *i m *2' one sees aroun( i m o = m i/2 the 
XX ~> ^ s— channel annihilation pole which would otherwise show up for higher values of tan/3 > 55 |r?j]. It 
affects significantly the fluxes from the Earth, where capture and annihilation are not in equilibrium figure 
[L6|c. The higher higgsino fraction also enhances the elastic cross section u x - q , but the stronger annihilation 
reduces the number of models which are both interesting from the cosmological point of view and in terms 
of Sun muon fluxes (figure p"s|c). As before, muon fluxes from the Earth are much too weak to be detected. 
Increasing Qewsb leads to the opposite effect: higher values of and reduced cosmological favoured region. 

Top mass: the current experimental value is m t = 174.3 ±5.1 GeV j26|. Our default choice up to now 
has been to adopt the central value. Lowering mt implies that the top Yukawa coupling Yt is smaller so that 
the RG decreasing of mj Iu is less efficient and the region where EWSB does not occur is wider (see eq. ^j, 
figure |l^ and |l^). Nevertheless the mixed higgsino region is larger favouring annihilation and capture, but 
there are no additional models combining a large i>/(i flux with an acceptable relic density (figure |l7], jl8] and 
p)[ ). Increasing the top mass leads to the opposite effect, which is less interesting for our study. It should 
also be noticed that varying Qewsb and/or m t allows for heavier neutralino within the fixed Qh 2 interval. 



7 Conclusion 

In this paper, we have analysed, within the framework of the constrained Minimal Supersymmetric Standard 
Model (a.k.a. mSugra), the possible indirect detection of neutralino dark matter through the neutrinos 
produced by its annihilation in the Earth and the Sun. In particular, we have presented an original study of 
the relative weight of various annihilation channels crucial to indirect detection. Fixing the halo dark matter 
density to p x = 0.3 GeV/cm 3 and allowing for the widest range of cosmological relic density 0.03 < Q.h 2 < 0.3, 
we find that a neutrino signal from the centre of the Earth is beyond reach of the present and even future 
neutrino telescopes. 

For the larger Sun, we find that a neutrino signal may be found, but only in the large mo "focus point" 
region |l^, where the neutralino has a larger higgsino component which increases both neutralino-nucleon 
elastic scattering and neutralino annihilation cross-sections. Since in this region, both are related by s- and 
i-channel crossing (exactly above tt threshold and approximately below) , any parameter change reducing the 
relic density inevitably increases the neutrino flux. When conversion into muons is taken into account with a 
5 GeV threshold on the muon energy, the only models providing both a relevant (but rather low) relic density, 
and fluxes around the current indirect detection sensitivity (10 3 /x km" 2 yr" 1 ) are those with m x < mt and 
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A0=0 ;tan(p)=10 ;^>0 
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m0 (GeV) mO (GeV) 

Figure 13: Branching ratio of the neutralino annihilation in the (jtoo: W.1/2) plane (same as figure ^ for tan/3 = 10). 
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Figure 14: a) Cosmologically favoured models for differents values of tan/3 in the (mc^m^) plane; subsetof these 
models generating a fj, flux from the Sun larger than b) 10 2 km -2 yr -1 ; c) 10 3 km -2 yr _1 . 



A0=0 ; tan(p)=45 ; n >0 ; mtop=174.3 GeV ; QEWSB= 1/2^ m s „ p1 m, lop2 




500 1000 1500 2000 2500 3000 

m0 (GeV) 



500 1000 1500 2000 2500 3000 

m0 (GeV) 



500 1000 1500 2000 2500 3000 

mO (GeV) 



a) 



b) 



c) 



Figure 15: a) Relic density, \i fluxes from the Sun (b) and the Earth (c) in the (mo, m 1 / 2 ) plane for Qewsb = lym^m^ 
and m t = 174.3 GeV. 



XX W + W~ or XX ~ ± * ZZ as dominant annihilation processes. Relaxing the tight mSugra constraints 
in specific directions could favorise such channels jl^]. The next generation of neutrino telescopes (with a 
km 3 size) will be much more efficient to probe mSugra models, especially for m x > nit ■ If a neutrino signal 
is found in this framework, then neutralino dark matter should also be accessible to future direct detection 
experiments, and a chargino with m x + < 350 GeV should be found in accelerators. 

To conclude, a comparison with previous related work is in order. Earlier studies Q ^ of neutralino 
indirect detection were performed in the unconstrained MSSM, where jj, is a free parameter, and mixed 
neutralinos occur rather easily. Taking an mSugra slice in this huge parameter space obviously introduces 
many correlations, for instance between the neutrino flux and the relic density or a\^!™ (and even o x c ^ l p ). 
The effect of this slicing in the {m x , u x c -p) plane was shown in [^] for small tan (3 and in for large tan (3. 
Our figure [] agrees with this last result, with the addition of the focus point region forming the upper-right 
cloud. We should note however that a comparison in the (mo, m^) plane is more difficult, specially at large 
tan/3, because of the theoretical uncertainties in the RGE used to translate these parameters into physical 
quantities. As an example, for a given value of tan /3, the slope of the EWSB boundary decreases between 
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Figure 16: a) Relic density, fi fluxes from the Sun (b) and the Earth (c) in the (mo, wi]^) plane for Qewsb = ^^/m^m^ 
and m t = 174.3 GeV. 
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Figure 17: a) Relic density, /i fluxes from the Sun (b) and the Earth (c) in the (mo,rn 1 /2) plane for Qewsb = ^/m^mj 
and m t = 170 GeV. 
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Figure 18: a) Relic density, \i fluxes from the Sun (b) and the Earth (c) in the (mcm^) plane for Qewsb = \^/mi 1 mi 2 
for m t = 170 GeV. 



pfifl , this work and fed]. As another example, the strong annihilation via s-channel A, Ho scalars appears 
above tan (3 > 50 in [ )5*l| and |5^ ], tan/3 > 35 in [Q, and tan/3 > 60 (or tan/3 = 45, Qewsb = ^/m^m^/5, 
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A0=0 ; tan(P)=45 ; u >0 




Figure 19: a) Dark matter favoured models for differents values of Qewsb and top mass in the (mo,m 1 /2) plane (a). 
Colors and symbols correspond to different values of mt and Qewsb'- e.g. (170 ; 1/2) means mt — 170 GeV and 
Qewsb = l^/m^m^. The models that moreover yield a muon flux from the Sun larger than > 10 2 km -2 yr -1 and 
> 10 3 km -2 yr 1 are shown in b) and c). 



see figure [16|) in the present work using SUSPECT. Finally, we find a muon event rate compatible with ]29[ , 
and larger than pl| . Part of this difference comes from a higher threshold (25 GeV instead of our 5 GeV). 
The rest might be attributed to a high sensitivity in the renormalization group equations at large mo and 
tan/3. It seems |J] that using SOFTSUSY §| or SUSPECT @ ^ may be safer than ISASUGRA§| in 
this region. 
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